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Abstract 

We consider for the time-dependent Maxwell's equations the inverse problem of 
identifying locations and certain properties of small electromagnetic inhomogeneities 
in a homogeneous background medium from dynamic measurements of the tangential 
component of the magnetic field on the boundary ( or a part of the boundary) of a 
domain. 
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1 Introduction 

Let be a bounded C 2 -domain in R , d = 2,3. Assume that contains a finite 
number of inhomogeneities, each of the form zj + aBj, where Bj C M. d is a bounded, 
smooth domain containing the origin. The total collection of inhomogeneities is B a = 
\JjLi{zj + aBj). The points Zj £ fi, j = 1, . . . , m, which determine the location of the 
inhomogeneities, are assumed to satisfy the following inequalities: 

\zj — Zji\ > co > 0, V j ^ f and dist(zj, dQ) > Co > 0, V j. (1) 

Assume that a > 0, the common order of magnitude of the diameters of the inhomo- 
geneities, is sufficiently small, that these inhomogeneities are disjoint, and that their 
distance to R d \ Q is larger than cq/2. Let an d £o denote the permeability and 
the permittivity of the background medium, and assume that /j,q > and eq > are 
positive constants. Let fij > and £j > denote the permeability and the permittivity 
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of the j-th inhomogeneity, Zj + aBj, these are also assumed to be positive constants. 
Introduce the piecewise-constant electric permittivity 

e a {x) = { £ ^ X ^ Q \\ . ( 2 ) 
\ £j, x G zj + aijj, j = 1 . . . m. 

If we allow the degenerate case a = 0, then the function £o(x) equals the constant £o- 
Assume that, the magnetic permeability is given by 

fJ-a(x) = /j-q, for all x £ fl. (3) 

Let v = v{x) denote the outward unit normal vector to f2 at a point on d£l, and dt = Jr. 

In this paper, we will denote by bold letters the functional spaces for the vector 
fields. Thus H s (tl) denotes the usual Sobolev space on f2 and H s (f2) denotes (H s (£l)) d 
and L 2 (f2) denotes (L 2 (Q)) d . As usual for Maxwell equations, we need spaces of fields 
with square integrable curls: 

H( curl ;Q)={ue L 2 (ft), curl u G L 2 (ft)}, 

and with square integrable divergences 

H(div ; n) = {u G L 2 (ft), div u G L 2 (ft)}. 

We will also need the following functional spaces: 

Y{9) = {u G L 2 (Q), div u = in Q}, X(Q) = H x (0) n y(fi), 

and TL 2 (<9f2) the space of vector fields on <9f2 that lie in L 2 (<9$7). Finally, the "minimal" 
choice for the magnetic variational space would be 

X N (n) = {v €H(curl;ft)nH(div;n); vxv = ondtt}. 

We use (•, •) for the duality bracket and (•, •) for the L 2 product. 

In the non flawed region Q an electric field/magnetic field pair (E,H) satisfies: 

curl E = -fJ, d t H infix (0, T), 

curl H = e d t E in n x (0, T), U 

— 1/2 

Let TH~rJ*[dU) denote the space of tangential vector fields on dQ, that lie in 
H-^idCl). The most common boundary data is 

(Exu)\ mmT) given in TH^ /2 (dfl). (5) 

By dividing the second equation in (|4|) by sq and taking the curl, we obtain in terms 
of the magnetic field: 

curl (— curl H) + fi d 2 H = in fi x (0, T), (6) 
and the boundary data is supposed to be given by 

(H x ^)Ux(o,T) = / given in TH^ /2 (dn). (7) 
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Moreover, we set 

H\ t=0 = if, d t H\ t=0 = ip in 0. (8) 

Here T > is a final observation time and tp,ip € C°°(Tl) and / G C oo (0,T;C oo (0fi)) 
are subject to the compatibility conditions 

dff\ t=0 = (AV) x Hen and c>f +1 /|t=o = (AV) x i/|an, 1 = 1,2,... 

Let i? a 6 M d be the magnetic field corresponding to the case of the presence of a 
finite number of small electromagnetic inhomogeneities. This field (under the assump- 
tion (J3j) ) satisfies 

' curl cur i H a ) + fi dfH a = infix(0,T), 
day (fx H a ) = infix(0,T), 

(9) 

H a \t=o = <P, d t H a \ t=0 = ip in O, 

. #a X ^|an x ( ,T) = /> 

It is well known that ([6]) has a unique solution H £ C°°([0, T]xO). It is also known (see 
for example [15]) that since Q is smooth (C 2 — regularity would be sufficient) the non ho- 
mogeneous Maxwell's equations (J9|) have a unique weak solution H a E C°(0, T; X(Q))n 
C 1 (0,T;L 2 (O)). Indeed, curl ff Q belongs to C°(0, T; n C 1 (0, T; L 2 (0)). 

Having found H a , we then obtain the field E a through the formula: 

dtE a = — curl H a . 

Our main goal in this paper is to determine, most effectively, properties of the in- 
homogeneities Zj + aBj, from over determined boundary information about specific 
solutions to ©. In particular, we study media that consist of a homogeneous(constant 
coefficient) electromagnetic material with a finite number of small inhomogeneities, 
and as our main result we derive asymptotic formulas for the perturbations in the 
(tangential) boundary magnetic fields caused by the presence of these inhomogeneities. 
Our formulas may be used to determine properties (location, relative size) of the small 
inhomogeneities in case a single, or a few (tangential) boundary electric fields and 
their corresponding(tangential)boundary magnetic fields are known. For stationary 
Maxwell's equations it has been known that the Dirichlet to Neumann map uniquely 
determines (smooth) isotropic electromagnetic parameters, see [14], |16j . [18] . We will 
provide in this paper a rigorous derivation of the inverse Fourier transform of a linear 
combination of derivatives of point masses, located at the positions Zj of the inhomo- 
geneities, as the leading order term of an appropriate averaging of (partial) dynamic 
boundary measurements of the tangential components of magnetic fields on part of the 
boundary. We refer the reader to [E] , [H] , [6] , and [8 J for discussions on closely related 
(stationary) identification problems. 

Our approach, aimed at determining specific internal features of an object based 
on electromagnetic boundary measurements, differs from pQ, [2], [3], [1], [20] and it can 
be regarded as constructive method, but until tested. 
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2 Asymptotic behavior 

We start the derivation of the asymptotic formula for curl H a x v with the following 
estimate. 

Lemma 2.1 The following estimate as a — > holds: 

\\dt{H a - -£0||l°c(o,T;L 2 (^)) + ll-^a ~~ H\\ L oa( 0;T . XN (n)) < Ca d , (10) 

where the constant C is independent of a and the set of points {z^JL-^ provided that 
assumption |7P holds. 

Proof. From ©-([9]), it is obvious that H a — H £ Xjv(f2), then due to the Green 
formula we have for any v € Xn(£1): 

[ fi d?(H a -H)-v+ [ — curl (H a —H)- curl v = V(— -— ) f curl H- curl v. 

(11) 

Let v a be defined by 

{ i ( 12 ) 

[ curl j- curl v a = dt{H a — H) in S7. 

Then, 

f — curl (fl- Q - H) • curl v a = - [ d t (H a - H) ■ (H a - H) = -\d t I \H a - H\ 2 
Jn £ a Jn * Jn 

and by Green formula, relation (|12[) gives: 

/ d 2 (H a — H) • v Q = J curl — curl <9fV Q ■ v„ 
Jn Jn £ a 

= — I — curl <9fV Q • curl v Q 
Jn £ a 

1 f 1 
= --#t / — I curl v Q | 2 . 

Thus, it follows from (1111) that 



/^o^t / — - 1 curl v a | 2 + dt [ \H a — H\ 2 = — 2Y^(— — ) / curl H ■ curlv a . 

Jn £ a Jn ~[ £ £ j Jzj+aBj 



Next, 



m 1 If 

I ^( ) / curl if- curl v«| < (7| | curl F|| L 2( Bo( ) || curl v a | Ils^). 

, i "» e J •'--.•.,/*, 

Since H G C°°([0,T] x O) we have 

m 

II curl#|| L2(Bci) < || curl H\\ Loo{Ba) a d Q2\ B j\)^ ^ C ° d ' 
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which gives 



1 1 



E( 

pi £ o £j J Zj + 



haBi 



curl H ■ curlv a | < Ca d \ \ curl v Q | |l2(h) 



and so, 



(14) 



fi d t [ — | curlv a | 2 +d t f \H a -H\ 2 <Ca d { [ — | curl v Q | 2 + f l^-H] 2 ) 1 / 2 . (13) 
Jn e a Jn isi £ a Jn 

From the Gronwall Lemma it follows that 

(/ l|curlv Q | 2 ) 1 / 2 + (/' \H a -H\ 2 )f/ 2 <Ca d . 
Jn £ a Jn 

Combining (|14p with the fact that 

\\dt(H a - -H") 1 1x00(0^^-1(0)) < C\\ curl v q ,|| LO o( 0)T;L 2(q)), 
the following estimate holds 

\\H a - ^||L°o(0,T;L 2 (n)) + \ \dt(H a - -ff)||L°°(0,T;L 2 (n)) < C« d - 

Now, taking (formally) v = dt{H a — H) in (jlip we arrive at 



(15) 



fi d t [ \\d t (H a -H)\ 2 +—\ curl (H a -H)\ 2 } =2V(---) / curl curl«9 t (tf 

Jn 1 £ a J -j £0 £j Jzj+ctBj 



By using the regularity of -ff in Q and estimate (|15p given above, we see that 
£(---) / ™rl curl d t (H a -H)\ < C\\ curl || H 2 (B J |<9t(i? a -tf) || H -i (P.) < 



3=1 

where C is independent of t and a, and so, we obtain 

1 



dt 



d t (H a - H)\ z + —I curl (H a - H) 



< Ca 



2d 



which yields the following estimate 

\\dt(H a - -£0||l°c(o,T;L 2 (^)) + W H <* ~ H\\L°°(o,T;X N (n)) < Ca d , 
where C is independent of a and the points {zj}'JL 1 . 

Now, we can estimate curl H a — curl H as follows. 

Proposition 2.1 Let H a and H be solutions to the problems (OJ) and (0) respectively. 
There exist constants < ao, C suc/i that for < a < oq the following estimate holds: 



□ 



curl (H a - ff)|Uoo ( o !T;L 2(n)) < Ca d , 



(16) 



Proof. To prove estimate (|16p it is useful to introduce the following function 



r T 

v(x) = / v(x,t)z(t) dt G L 2 (n), 
Jo 

where v G L 1 (0, T; L 2 (£l)) and z(t) is a given function in Cq°(]0,T[). 
Then, 

H(x) = [ H(x,t)z(t) dt and H a (x) = [ H a (x, t)z(t) dt G X(D,), 



(17) 



which by relation (jlOp give 

f (# a -#) EH 1 )!!), 

curl curl (H a - H ) = 0(a d ) in 0, 
div (H Q - J?) = in 0, 
k (H a - H) x z^|an = 0, 

and so, 

|| curl | | L3(n) =0(a d ). 

The fact that curl (H a - H) belongs to L°°(0, T; L 2 ($7)) and by using (TTT]) and 
we arrive at: 

/ | curlF a (x,t) - curl H(x,t)\ 2 dx = 0(a 2d ) a.e. in t G (0,T), 
in 



(18) 



which means that 



curl (H a - JET)|| L a(n) = 0(a d ) a.e. in t G (0,T). 



This equation can be bounded easily according to t G (0, T). Thus, estimate ([To]) holds. 



□ 



Before formulating our main result in this section, let us denote j = 1, . . . ,m 
the unique vector-valued solution of the following free space Laplace equation: 

' A$j = in Bj, and ~R d \ B~, 
$j is continuous across dBj, 



£0 



<9za- 



(19) 



lim |$-(y)|=0, 



»+oo 



where denotes the outward unit normal to dBj, and superscripts — and + indicate 
the limiting values as the point approaches dBj from outside Bj, and from inside Bj, 
respectively. The existence and uniqueness of this <&j can be established using single 
layer potentials with suitably chosen densities, see [6] for the case of conductivity 
problem. For each inhomogeneity Zj + aBj we introduce the polarizability tensor 
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Mj which is a d x d, symmetric, positive definite matrix associated with the j-th 
inhomogeneity, given by 

(M,-) M = e k • ( / (uj + - l)^\+(y))y ■ ei daj(y)). (20) 

Here (ei, . . . , e^) is an orthonormal basis of M. d . In terms of this function we are able 
to prove the following result about the asymptotic behavior of curl H a ■ Vj\d( Zj +aBj)+- 

Theorem 2.1 Suppose that (QJ) is satisfied and let <&j,j = l,...,mbe given as in \19\). 
Then, for the solutions H a , H of problems |5j) and respectively, and for y £ dBj 
we have 

( curl H a (zj + ay) ■ Vj)\a( Zj +aB j )+ = curl H(zj,t) ■ Uj (21) 
+ (1-^)|^|+(V)- curlF(^)+ (l). 

£0 OVj 

The term o(l) uniform in y £ dBj and t £ (0, T) and depends on the shape of {Bj} 1 ^-^ 
and Q, the constants cq, T, £q, {£j}JLi, the data <p,ip, and f, but is otherwise inde- 
pendent of the points {zj}™^. 

Proof. 

Let TC a = curl H a (x, t) and TLq = curl H(x, t). Then, according to ©-(UJ) we have 

HodfH a — curl — H a = and curl H a = 0, for (22) 

We restrict, for simplicity, our attention to the case of a single inhomogeneity, i.e., 
the case m = 1. The proof for any fixed number m of well separated inhomogeneities 
follows by iteration of the argument that we will present for the case m = 1. In order 
to further simplify notation, we assume that the single inhomogeneity has the form 
aB, that is, we assume it is centered at the origin. We denote the electromagnetic 
permeability inside aB by and define the same as <I>j, defined in (|19p . but with 
Bj and £j replaced by B and £*, respectively. Define v to be the outward unit normal 

to dB. Now, following a common practice in multiscale expansions we introduce the 

x SI 
local variable y = — , then the domain (l = (— ) is well defined. 

a a 
Next, let w be given in C(J°(]0, T[). For any function v S L 1 (0, T; L 2 (f2)), we define 

f T 

v(x) = / v(x, t) zu(t) dt £ L 2 (J7). 
Jo 

We remark that dtv(x) 
satisfies 

< 



T 

v(x, t)zu'(t) dt. So that we deduce from (f2"2"|) that H c 



o 



1 " f 
curl — TC a = / H a w (t) dt in O, 

jo 

curl H a = in fi. 
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Analogously, H satisfies 



to 



curl TL 



H w (t) dt in ft, 



o 



curl H = in ft. 



Indeed, we have H a x u = H x v = curl qq/ x v on the boundary <9ft, where curl 
is the tangential curl. Following [3] and p], we introduce q* a as the unique solution to 
the following problem 

f Q _ 

Aq* a = in ft = (-) \ B and in B, 

a 

q* a is continuous across dB, 

dq* dq* 
e —^\ + - = ~(eo - e*)H(ay) ■ v on dB, 

av ov 

q* = on <9ft. 



The jump condition 



0q 



£0-^1+ - = ~(eo - £*)H(ay) ■ v on 95 



<9z/ 



<9z/ 



guarantees that TL a (x) — TL(x) — grad ^Q^C-) belongs to the functional space Xn(£1), 
where grad qq is the tangential gradient. Since 



I ^ * x 
curl — (H a -H- grad q* a (-)) 



a: 



# Q - x(0 \ a#)# + — x(atB)H] w"{t) dt in ft, 



curl (7Y a — H — grad — )) =0 in ft 



a 



k (H a — H — grad Xi/ = on 9ft, 

where is the characteristic function of the domain uj, we arrive, as a consequence 
of the energy estimate given by Lemma l2.lt at the following 

f (H a -TL- grad y q* a (§)) € X N (Sl), 

L x 
curl — (Ha — H— curl y q* a ( — )) = 0(a) in f2, 

^ x 
curl (H a — TL — grad v ql(— )) =0 in ft, 
y a 

I (H a - W - grad „<£(*)) x i/ = on aft. 
From [4] we know that this yields the following estimate 



curl — (TL a - H- grad y g* f 



a 



| L 2 ( n) + ||H Q — H— grad ^g* (-)|| L 2 (C) < Ca, 



a 



and so, 



(H a - H - grad <?* (-)) • z/| 



0(a) ond(a5). 
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Now, we denote by q* be the unique (scalar) solution to 
'Ag*=0 mM. d \BandmB, 



is continuous across dB, 



dq* . dg* 

lim g* = 0. 

\y\-*+oo 



-(e - e m )H(0) ■ v ondB, 



In the spirit of Theorem 1 in [6] it follows that 

||( grad y q* - grad y q*)(~)\\jp(n) < Ca 1/2 , 

which yields 

(Wo — Ji— grad «<?*(—)) ■ ^ = o(l) on d(aB). 



a 



Writing q* in terms of $>* gives 
T 

L 



curl H a (ay)-v)\ d(aB) +-u- curl H(0, t)-{— curl H(0,f) (it = o(l), 



for any ot G Cq°(]0, T[), and so, by iterating the same argument for the case of m 
(well separated) inhomogeneities Zj + aBj,j = l,...,m, we arrive at the promised 
asymptotic formula (f2Tj) . 



□ 



3 Reconstruction 

Before describing our identification and reconstruction procedure, let us introduce the 
following cutoff function (3{x) £ Cg°(f2) such that (3 = 1 in a subdomain D,' of f2 that 
contains the inhomogeneities B a and let rj E We will take in what follows H (x, t) = 
where j? -1 " is a unit vector that is orthogonal to r] which corresponds to 
taking <p(x) = rj ± e ir, ' x , ip(x) = —i^/eo\n\r]- L e vri ' x , and f(x,t) = n 1 - x ^"^^v^oN* and 
assume that we are in possession of the measurements of: 

curl H a xu on T x (0,T), 

where T is an open part of <9f2. Suppose now that T and the part V of the boundary d£l 
are such that they geometrically control 0, which means that they satisfy the geometric 
control hypothesis of the work of Bardos, Lebeau and Rauch in [5]: 

Definition 3.1 Let T be an open subset of dQ and T a positive number. One says 
that (r,T) geometrically control £1 if for every geometrical optic ray s ^ 7(s), there 
exists so e]0,T[ such that j(sq) € Tx]0,T[ and 7(so) non diffractive point. 

It follows from [15] (see also [11], [9] and [10]) that we can construct (a unique) 
g v E H^(0,T;TL 2 (T)) (by the Hilbert Uniqueness Method) such that the unique weak 
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solution w v to 



' (<9 t 2 + curl curl )w r) = in fi x (0, T), 
divw,, = infix (0,T), 
^| t=0 = (3(x)n ± e ir i' x , d t w v \t=o = in fi, 



(23) 



^ x ^bc\rx(o,T) 



0, 



k x i/|r X (o,T) = 

satisfies w n {T) = d t w v (T) = in fi. 

Let 9 V e ^(O^TL 2 ^)) denote the unique solution of the Volterra equation of 
second kind 



d t e ri {x,t) + 



-i\v\(s-t) 



9 v (x,0)=0 forxer. 



(9 v (x, s) — i\r]\dt0 v (x, s)) ds = g v (x, t) for x £ T, t G (0, T), 



(24) 



The existence and uniqueness of this 9 V in H 1 (0, T; TL 2 (T)) for any rj 6 R rf can be 
established using the resolvent kernel. However, observing from differentiation of 
with respect to t that 6^ is the unique solution of the ODE: 



d 2 9 v -0 V = eWdtie-Wgr,) for x E T, t 6 (0, T), 
e,,(x,0) = 0,d t 9 v (x,T) = forxGT, 



(25) 



the function ^ may be find (in practice) explicitly with variation of parameters and it 
also immediately follows from this observation that 6 V belongs to H 2 (0, T; TL 2 (T)). 
We introduce v v as the unique weak solution (obtained by transposition as done in 
[13] and in p] [Theorem 4.2, page 46] for the scalar function) in C°(0,T;X(fi)) n 
C^TjL^fi)) to the following problem 

' (<9 t 2 + curl curl )v v = in fi x (0,T), 
div^ = infix (0,T), 
^=0 = in fi, 



d(zj+aBj) 



B ,)€Y(Q) infi, 



9t^| t=0 = - — )?? x (i/j + (— 

^ X ^bnx(0,T) = 0. 
Then, the following holds. 



Proposition 3.1 Suppose that T and T geometrically control fi. For any rj 6 M. d and 
r] 1 - unit vector in M. d that is orthogonal to r\, we have 



/ / g v ■ ( curl v v x is) = a d y~^ £ (1 
Jo Jt 



£0 



rj x 



(26) 



>>>)</</) •/r L ^,(//) + 0(n"). 
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Proof. Multiply the equation (<9 2 + curl curl )v v = by w v and integrating by parts 
in t G (0,T), we get 

/ / i^t + cur ^ cur l ) v v w n = / cup l cur ^ v v w v + / / dt v v w ri 
Jo J n Jo Jn Jn Jo 

/ curl curl t),w, + / %^™?7lt=0 - dtV^I*^ _ / / dtv v d t w v 
o ./n in Jn Jo 

= / curl curlv„w,,+ / d t v v w v \ t =o + / ^^^^=0+ / / v v d 2 w v . 
Jo Jn Jn Jn Jn Jo 

So, by Green's formula, 

[ T [ (d 2 + curl curl )v v w v = -a d ~ x - ^-)e 2 ^^r ? x ( / ( Vj + 

Jo Jn j=i £ ° J dB -> 



[ £ J- - l)^l\ + (y))e^'y) ■ ^r, 1 ds 3 {y) 
-Bq 1 J j g n • ( curl x v) = 0. 



Therefore 



a d-i Vi(l - ^)e 2 ^r ? x ( / + (fi - i)_^| + ( y )) e ^) • /3( y )^ d Sj (y) = 

~ e o 1 / / 9v ' ( curl ^ x v )- 
Jo Jr 

Now, we take the Taylor expansion of a d ~ 1 e iar, ' y in the left side of the last equation and 
we use the definition of the cutoff function f3(x), we obtain the convenient asymptotic 
formula (|26j) . 

□ 

To identify the locations and certain properties of the small inhomogeneities B a let 
us view the averaging of the boundary measurements 

curl H a x Hrx(o,T), 

using the solution 9 n to the Volterra equation or equivalently the ODE ([25]) , as a 
function of rj. The following holds. 



Theorem 3.1 Let rj 6 M. d and n 1 - be a unit vector in M. d that is orthogonal to n. Let H a 
be the unique solution inC°(0,T;X(n))<lC 1 (0,T;L 2 (n)) to the Maxwell's equations 
with <p(x) = r/- L e jr? ' x , tjj(x) = — i^fe^r}\r)^ e %ri ' x , and f(x,t) = r] ± e l ' n ' x ~' l ^^ t . Suppose 
that r and T geometrically control Q, then we have 

J J ' ( cur l H a x v — curl H x v) + dtOr) ■ dt{ curl H a X v — curl H x v) 

m 

a d Y,(Bo ~ B^e^iri x Mj(n)) • ^ + 0(a d ), 

(27) 
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where 9„ is the unique solution to the Volterra equation A25\) with g„ defined as the 
boundary control in f23\) and Mj is the polarization tensor of Bj, defined by 



(28) 



JdBj £ ° v j 



Here (ej, • • • , e^) is an orthonormal basis ofM. d . The term 0(a d ) is independent of the 
points {zj, j = 1, • • ■ , to}. 



Proof. From dtO„(T) = and ( curl H a xv — curl H x f)|(=o = the term 
dt{ curl H a x v — curl H x v) has to be interpreted as follows 




o Jr 



d t e v - 



ndtO v -dt{ curl H a x v— curl Hx v) = — / / (9 t 2 ^ • ( curl if Q x v — curl # x v). 
Jo Jr 

(29) 

Next, introduce 

H a , v (x, t) = H(x, t)+ f e-^iK^x, t-s)ds,xen,te (0, T). (30) 
Jo 



We have 

cT 




>o Jr 

fT 



6 V ■ ( curl H a x v - curl H x v) + • <9 t ( curl x v — curl x 



o Jr 



fT 



+ 



JT 



% t ■ ( curl H a xv— curl fl^Q, fl X v) + <9t^ • <9t( curl H a x u — curl H a>T) x u) 



'^"Knix, t-s)xuds + d t 9 v ■ d t I e'^^v^x, t - s) x v ds 



Since 8 V satisfies the Volterra equation (|2"5|) and 

d t { P e- l ^ Ms v v (x,t - s)xu ds) = dtie-^W f e^^n^x, s) x v ds) 



o 



J v (x, s) x v ds + 1*77(2;, t) x u, 



we obtain by integrating by parts over (0, T) that 



t-T 




Jr 



-V^H^( X) t _ s ) x v ds + d t 6 v d t e 
(v v (x, t)xv)- (8 t 9 v + / ^(aJeV/^K*--) ds) 



^(x, t — s) x 1/ ds 



-i-s/^Hie-^^dte^t)) ■ / e^'I'^Cs.a) x v ds dt 



rT 



■j v ( x ,t) x v ■ (d t e v + J (e v (s) - i^UdtO^y^^ ds) 



dt 



Jr 



g v (x,t) ■ ( curl v v (x,t) x v) dt 
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and so, from Proposition 13.11 we obtain 



curl H a x v — curl H x u) + dtO v ■ dt( curl H a X u — curl H x v) 



a 



+ 



E(l - % 2 "^ (17 X ( / fa + - l)^|+(y))r? • »)) • ^ cfyfo) 




./r 



^ • ( curl H a xv— curl H a>r) x v) + 3^ • <9i( curl ff a 



x v 



— curl Hq, )7? x v) + o(a ). 
In order to prove Theorem 13, II it suffices then to show that 

Jo Jr L 



6 v -( curl H a xv— curl H a>Tl xv)+dt6 ri -dt( curl H a xv— curl H a ^xv) = o(a d ). 

(31) 



Since 



( 1 /■* 

(3 t 2 - curl — curl )(/ e'^^v^x^ 
£ o Jo 



s) ds) 



£<(1 - x ("i + (- - l)7T^I + (y)) ei,? ' 2: '^+ aB3 )e-^ |,?l * m fi x (0,T), 



( / e -^N^( X) t_ s ) ds )| t=0 = o,9 t ( 



( / e -VSShl^( S) t _ s ) dfl ) x u\ dnx{0)T) = 0, 



s v v (x, t — s) ds)\t=o = in f2, 



< Jo 

it follows from Theorem 12.11 that 

(d t 2 - curl — curl )(H a - H a „) = 0(a d ) € Y(Q), in ft x (0, T), 

£0 

(H a — H ajV )\ t= o = 0,dt(H a — H a ^)\ t= Q = in Q, 



y {H a — H a>r) ) x v\ 



dnx(0,T) 



0. 



Following the proof of Proposition 12.1^ we immediately obtain 

\\H a -H a J\ L 2 (n) = 0(a d ),te (0,T),x G O, 

where 0(a d ), d = 2,3 is independent of the points {zj}™^- To prove (|3T|) it suffices 
then from (|29p to show that the following estimate holds 

|| curl H a x u - curl H a ^ x v\ \ l ^(o,t-,tl^{T)) = 0(a d ). 
Let 9 be given in Cg°(]0,T[) and define 



H a , v (x)= [ H a , n {x,t)6{t) 
Jo 



dt 
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and 

H a (x) 

From definition (1301) we can write 



H a (x,t)9(t) dt. 



{H a — H a ) G H 1 (17), 

curl curl (H a - H a ) = 0(a d ) G F(O) in fl, 
div (fl" a - fl" a ) = in Q, 
(H a - H a ) x v\dn = 0. 



(32) 



In the spirit of the standard elliptic regularity [7] we deduce for the boundary value 
problem (f32l) that 

|| curl (H a - H a ) x z/|| L 2 (r ) = 0(a d ), 
for all 9 G Cg°(]0,T[); whence 

|| curl (H a - H a ) x ^|| L 2(r) = 0(a d ) a. e. in t G (0, T), 

and so, the desired estimate (|27|) holds. The proof of Theorem 13.11 is then over. □ 

Our identification procedure is deeply based on Theorem 13.11 Let us neglect the 
asymptotically small remainder in the asymptotic formula (|27|) . and define ^aiv) by 



M»7) 



r 



0„ • ( curl -ff)x!/) + • d t { curl (ff Q - H) x u) 



Recall that the function e 2ir, ' Zj is exactly the Fourier Transform (up to a multiplicative 
constant) of the Dirac function e)_2 2 (a point mass located at — 2zj). From Theorem 
13.11 it follows that the function e 2ir i' z i is (approximately) the Fourier Transform of a 
linear combination of derivatives of point masses, or 



3=1 



2z, 



where Lj is a second order constant coefficient, differential operator whose coefficients 
depend on the polarization tensor Mj defined by (f28|) (see [6] for its properties) and 
H a (r/) represents the inverse Fourier Transform of H a (r]). 

The method of reconstruction consists in sampling values of ^a(v) a ^ some discrete 
set of points and then calculating the corresponding discrete inverse Fourier Transform. 
After a rescaling the support of this discrete inverse Fourier Transform yields the loca- 
tion of the small inhomogeneities B a . Once the locations are known we may calculate 
the polarization tensors (Mj) 1 J^ =1 by solving an appropriate linear system arising from 
(1271) . This procedure generalizes the approach developed in [3] for the two-dimensional 
(time-independent) inverse conductivity problem and generalize the results in [1] to 
the full time-dependent Maxwell's equations. 
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